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Abstract 

The existence question for two-dimensional symmetric steady waves travelling 
on the surface of a deep ocean beneath a heavy elastic membrane is analyzed as a 
problem in bifurcation theory. The behaviour of the two-dimensional cross-section 
of the membrane is modelled as a thin (unshearable) , heavy, hyperelastic Cosserat 
rod, and the fluid beneath is supposed to be in steady two-dimensional irrotational 
motion under gravity. When the wavelength has been normalized to be 2n, and 
assuming that gravity and the density of the undeformed membrane are prescribed, 
there are two free parameters in the problem: the speed of the wave and drift 
velocity of the membrane. 

It is observed that the problem, when linearized about uniform horizontal flow, 
has at most two independent solutions for any values of the parameters. When 
the linearized problem has only one normalized solution, it is shown that the full 
nonlinear problem has a sheet of solutions comprised of a family of curves bifurcating 
from simple eigenvalues. Here one of the problem's parameters is used to index a 
family of bifurcation problems in which the other is the bifurcation parameter. 

When the linearized problem has two solutions, with wave numbers k and I such 
that max{fc, I}/ min{fc, 1} ^ Z, it is shown that there are three two-dimensional 
sheets of bifurcating solutions. One consists of "special" solutions with minimal 
period 2n/k; another consists of "special" solutions with minimal period 2n/l; and 
the third, apart from those on the curves where it intersects the "special" sheets, 
consists of "general" solutions with minimal period 2ir. 

The two sheets of "special" solutions are rather similar to those that occur 
when the linearized problem has only one solution. However, points where the 
first sheet or the second sheet intersects the third sheet are period-multiplying (or 
symmetry-breaking) secondary bifurcation points on primary branches of "special" 
solutions. This phenomenon is analogous to that of Wilton ripples, which arises in 
the classical water-wave problem when the surface tension has special values. In the 
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case of Wilton ripples, the coefficient of surface tension and the wave speed are the 
problem's two parameters. In the present context, there are two speed parameters, 
meaning that the membrane elasticity does not need to be highly specified for this 
symmetry-breaking phenomenon to occur. 

Keywords: hydrodynamic waves, hydroelastic waves, nonlinear elasticity, free boundary prob- 
lems, travelling waves, bifurcation theory, secondary bifurcations, Wilton ripples, Lyapunov- 
Schmidt reduction, symmetry-breaking. 
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1 Introduction 

The existence question for symmetric, 27r-periodic steady waves travelling with speed 
Co on the surface of a heavy, inviscid fluid which is at rest at infinite depth beneath a 
heavy, thin (unshearable) elastic membrane was considered in [7] as a global problem in 
the calculus of variations. Here we use local methods to study the bifurcation of such 
waves. We will show that there are two free parameters, cq (the wave speed) and d 
(the membrane drift velocity) and that, when the problem is linearized about a uniform 
horizontal stream with the membrane unstretched, there are no more than two linearly 
independent solutions. (See equation (|5.2|) , in which Ai = p(co — d) 2 , A2 = Cq, p is the 
density of the membrane and C is the Hilbert transform of a 27r-periodic function.) When 
there is no non-zero solution, nonlinear waves do not bifurcate from uniform horizontal 
streams; when there is only one solution, a sheet of solutions representing a parameterized 
family of bifurcations from simple eigenvalues occurs; when there are two independent 
solutions, there bifurcate three sheets of small-amplitude periodic waves. The latter 
corresponds to the presence of secondary bifurcations from curves of "special" solutions, 
the hydroelastic analogue of what are known as Wilton ripples [9j, as described in the 
Abstract and in Section [L2"l To quote from 0], "Waves characterized by two dominant 
modes are often called Wilton's ripples in the literature in reference to Wilton's paper 
(1915). It turns out that the phenomenon described as Wilton's ripples was accounted 
for at least twice prior to Wilton's paper : in an unpublished addendum to an essay that 
Bohr (1906) wrote in order to win the Royal Danish Academy prize on the theme 'The 
surface tension of water' and in a paper by Harrison (1909)." (Bohr's essay is [3].) 

A key feature of the present analysis is the reduction of the physical problem (|1.3I a-g) to 
an equation (14. 7[) for one 27r-periodic function of one real variable and two parameters. 

Physical Problem 

This problem is described in detail in [7]. To summarise, we consider waves on the surface 
of an infinitely deep ocean beneath an elastic membrane under the assumption that there 
is no friction between the membrane and the fluid. Since the water depth is infinite, there 
is no loss (after normalizing length scales appropriately) in restricting attention to waves 
that have period 27r in the horizontal direction. The fluid's Eulerian velocity field is 
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supposed to be two-dimensional and stationary at the same time as the material of the 
membrane is in motion, driven by gravity, by forces and couples due to its elasticity and 
by pressure from the fluid. The resulting mechanical behaviour of the surface membrane 
is modelled by regarding its cross-section as a heavy, unshearable, hyperelastic Cosserat 
rod, using the treatment in Antman [TJ Ch. 4]. We deal with this first. 

Membrane Elasticity. Let (x, 0) G K 2 be the rest position of a material point in 
the membrane cross-section and let r(x) G M 2 be its position after deformation. In the 
notation of Antman, 'd(x) is the angle between the horizontal positive semi-axis and the 
vector r'[x) (where ' means d/dx). Let 

v(x) = \r'(x)\ and fi(x) = $'[x). (1.1) 

Thus v(x) is the stretch of the membrane at the point r(x) and 

(j,{x) 



a(r(x)) 



v(x) 



is its curvature. We suppose that the elastic properties of the membrane are described 
as follows. 

HI. (Hyperelasticity) There exists a C°°-function E(is,fj.) > 0, v > 0, fi G M, such 
that, after the deformation (a;, 0) i— > r(x), the elastic energy in the reference segment 
{(x, 0) : x G [2:1,2:2]} is 

£(r) = / E(v(x), fi(x)) dx, 
where v{x),ii{x) are defined in (jl.ip . E is called the stored energy function. □ 



We also assume that the reference configuration, unstretched and unbent, is a local 
minimum of the elastic energy, which is locally convex. 

H2. (Rest state and local convexity) 

^(1,0) = S 1 (1,0) = S 2 (1,0)=£; 12 (1,0) = 0, S n (l,0)>0, £ 22 (1,0)>0. 
Subscripts 1, 2 denote partial derivatives with respect to v, /i, respectively. □ 

Remark. Since this is a study of bifurcating waves, regularity questions that are sig- 
nificant for large-amplitude waves are unimportant here. It is therefore for convenience 
only that we suppose E G C°°. With a little more technical effort the theory can be 
developed for E with much less regularity. □ 



Travelling Waves. For a periodic travelling wave, the position at time t of the material 
point with Lagrangian coordinates (x, 0) in the undeformed membrane is assumed to be 
given by 

R(x, t) :— (x + dt + u(x — ct), v(x ~ ci)) , 
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where u and v are 27r-periodic and c, d £ R. Let Co = c + rf. Then the surface profile at 
time i is the curve 

<S t = {(x + dt + u(x - ct), v(x - ct)) :i£lj 
= {(s + u(s),v(s)) :sel} + M,0) 
=: 5+ (co*,0). 

Thus St is represented by a profile S of fixed shape propagating from left to right, say, 
without changing shape at a constant velocity cq while at the same time the material 
point with Lagrangian coordinates (x,0) has temporal period 2n/c relative to a frame 
moving with speed d. We refer to Co as the wave speed and to d as the drift velocity of 
the membrane, both calculated relative to the fluid at rest at infinite depth. 

Since the membrane is in motion relative to the moving frame, there are inertial effects 
due to its mass, but it is supposed throughout that there is no friction between the fluid 
and the membrane. Under this assumption, it was shown in [7j that the inertial effects 
lead to an equivalent steady- wave problem in which the wave speed and the drift velocity 
coincide, and the stored energy function is perturbed by a quadratic term. 

Remark. To motivate this observation and what follows, it is worth observing the corre- 
sponding situation that arises when travelling- wave solutions R to an analogous nonlinear 
wave equation 

(E'(R X )) X = pR tt , p>0, 

are sought. This equation describes longitudinal motion in a one-dimensional elastic rod 
for which E is the stored energy function, p is the undeformed density and R(x, t) G R is 
the position at time t of the point with Lagrangian coordinate x S R. It is significantly 
simpler than the hydroelastic wave problem because here there are no body forces and 
because the shape of the rod does not change, only the relative positions of its material 
points in a straight line change with time. 

First note that a stationary (time-independent) solution satisfies (E'(R X )) X = and is 
given by a critical point of the potential energy functional 

/ E(R x )dx. 
Jo 

On the other hand, R is a periodic travelling wave with drift velocity d if 

R(x, t) = cot + r(x — ct), 

for some c and cq, where r(s + 2n) = 2n + r(s) and d = cq — c. Here we regard the variable 
s ( = x — ct) as a steady Lagrangian coordinate for travelling waves. The equation for 
r(s) is then 

{E'(r s )-c 2 pr s } s = 0, (1.2) 
which corresponds to critical point of 
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Therefore periodic travelling waves correspond to a boundary-value problem for station- 
ary solutions of a nonlinear wave equation with a different stored energy function 

E(p) := E(p) - |cV, 

instead of E(p). In |7J, E_ is called the pseudo-potential energy of travelling waves. If c 
is large E_ is not convex, even when E is convex [6] . 

Note that if r and c correspond to a travelling wave, then a family, parametrized by 
Co G R, of travelling waves with drift velocity d, is given by 

R(x, t) = cot + r(x — ct), d = Co — c. 

However, in the hydroelastic wave problem, the interaction of the membrane with the 
fluid means that the dependence of waves on both parameters is not so trivial. In fact, 
we will see that the solutions r(s) = s, Co and c arbitrary, to this one-dimensional 
problem, when combined with a wave profile with zero elevation, is a family of trivial 
solutions of the hydroelastic wave problem which we describe next. They correspond to 
an undeformed membrane drifting with velocity d = Co — c on the surface of a uniform 
flow with horizontal velocity Co- □ 

To summarise the analogous treatment in [7] of the equivalent hydroelastic travelling- 
wave problem, let (s,0) be the steady Lagrangian coordinate, r(s) its deformed position 
and let p and v be defined as in (jl.ip . with s in place of x. Let p be the density of the 
undeformed membrane section and let rj(s) = j • r(s), where j is the unit vector in the 
upward vertical direction (77 is the wave elevation). In pj eqn. (1.8)] it is shown that v, 
p and rj satisfy 

— c 2 v{sf+p{s)E 2 {v{s), l x{s))-E{ V ,p)-gfy n {s)] = 0, (1.3a) 

which is the analogue of (|1.2|l in the present situation. The pressure P in the fluid, 
internal forces and gravity combine to deform the membrane. Thus, from [2 (1.7e)], 

p i r ) = - m) — c pv) H , (1.3b 

where, as in [7] , we assume that the material in one period of the membrane surface is a 
deformation of an interval of length 2-k of the reference membran^l: 

S n ([0, 2tt] Xl) = {r(s) : s E [s , s + 2tt}} for some s e R. (1.3c) 

Fluid Motion. In the moving frame the fluid velocity field is two dimensional, stationary 
and irrotational, the membrane cross-section coincides with a streamline and the flow at 
infinite depth is horizontal with velocity — Co- The surface S is a zero level line for the 

3 The periodic wave does not require additional membrane mass as it passes underneath. 
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stream function ip(X,Y), which is harmonic and represents horizontal laminar flow at 
infinite depth. So 

Aip = below S, (f .3d) 

ip — on S (the kinematic boundary condition), (f .3e) 

Vtp(X, Y) -> (0, c ) as Y -> -oo. (1.3f) 

Membrane-fluid interaction. The dynamic boundary condition takes the form 

-^\\7^(X,Y)\ 2 ~gY+^=P(r( S )) when (X,Y) = r(s) E S. (1.3g) 

Here P is given by () 1 .3b[) and the left side of Q1.3g[ ) is the pressure in the fluid. 

1.1 The Free-Boundary Problem 

A steady hydroelastic wave with speed Co and drift velocity d is a non-self-intersecting 
smooth 27r-periodic curve S in the plane for which there exists a solution of (jl.3p with 
c = Co — d. Since we are interested in symmetric waves, we require that S is symmetric 
about a vertical line. To examine the solution set of this elaborate system we will use 
standard bifurcation-theory methods based on the implicit function theorem. Before 
going into the details, we give a schematic outline of what we have achieved (much 
remains to be done) . Throughout we treat g (gravity) and p (density of the undeformed 
membrane) as constants and regard 

A i: =c 2 P , A 2 :=cg, A= (Ai,A a ) eR 2 , 

as the physical parameters. (Because bifurcation theory studies the existence of solutions 
with small amplitudes and slopes, the question of self-intersection of wave surfaces does 
not arise in the present study. By contrast, in the theory of large-amplitude waves [2], 
considerable effort is required to ensure that there is no self-intersection.) 

1.2 Bifurcation Picture 

In this section we explain schematically the geometric picture of small amplitude hy- 
droelastic waves close to a bifurcation point. See Theorems 17.11 and 18.11 for a detailed 
statement. 

The first observation is that, for all choices of the two independent parameters, the 
problem, when linearized at the trivial solution of uniform horizontal flow under an 
undeformed membrane, has at most two linearly independent solutions. If there is only 
one linearized solution when (Ai, A2) = (A£, A£) say, there is at most only one linearized 
solution for all nearby (Ai, A2). Therefore, with either one of the parameters held fixed, 
there are bifurcations from simple eigenvalues with respect to the other parameter. Their 
union is a two-dimensional sheet of solutions which bifurcates from (A*, A|). The details 
are in Section [3 
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On the other hand, suppose that at A* = (A*, A|) the linearized problem has two solutions 
cos(fcr) and cos(Zr), where k and I are positive integers with max{fc, I}/ min{fc, 1} 
Z. By restricting attention to solutions in Zj~ = span {cos(j kr) : j £ N}, or in Zi = 
span {cos(j7r) : j 6 N}, the problem may be reduced to one of "bifurcation from a simple 
eigenvalue" for particular solutions that have minimal period 2tt /k or 2tt/1, respectively. 
This is straightforward and similar to what is done in Section [7] Locally we obtain a 
sheet of solutions of minimal period 2n/k and a sheet of solutions of minimal period 
2n/l. Each of these sheets is locally the graph of a function which gives A2 in terms of 
the wave amplitude and Ai (the roles of \ 2 and Ai can be reversed). We will refer to 
these solutions with minimal periods less than 2ir as "special" solutions. 

In Section [5] we show that in addition to these two sheets of "special" solutions there is 
a two-dimensional sheet of "general" solutions and that the sheet of "general" solutions 
intersects each of the sheets of "special" solutions in a curve. The solutions on the sheet 
of "general" solutions, except where it intersects the sheet of "special" solutions, have 
minimal period 2tt. Therefore the general solutions on this sheet represent a symmetry- 
breaking (or period- multiplying) secondary bifurcation on the curves of special solutions. 
This is the hydroelastic analogue of Wilton ripples, a type of water wave which arises 
in the presence of surface tension [3 [9] . In Wilton-ripple theory there are also two 
parameters, the wave speed and the surface tension coefficient (which measures surface 
elasticity). Wilton ripples bifurcate from uniform streams at certain values of the wave 
speed, when the surface tension has particular values. Here the two parameters are 
independent of the elasticity of the membrane. Therefore the wave and drift speeds can 
conspire to produce ripples, for any prescribed elastic membrane. 

After Lyapunov-Schmidt reduction, when the linearized problem has two independent 
solutions, cos(kr) and cos(Zr), there are two bifurcation equations in four unknowns 

$ fc (<i,i 2 ,Ai,A 2 ) = 0, $i(*i>*2,Ai,A 2 ) = o, (1.4) 

where t\ and t 2 near are the coefficients of cos(fcr) and cos(ir), respectively. Solu- 
tions with t\ = or t 2 = correspond to "special" solutions. The hypothesis that 
max{fc, I}/ min{fc, 1} ^ Z leads to the key observation that, for all ti,t 2 near 0, 

$ fc (0, t 2 , Ai, A 2 ) = 0, ^i(t u 0, Ai, A 2 ) = 0. 

Then each of the sheets of "special" solutions is found by solving one equations in three 
unknowns, as is done in Section [7l 

To find the "general" solutions, we seek solutions of (| 1 .4[) with neither t\ nor t 2 equal to 
0. This problem is reduced in Section [8] to a desingularized one of the form 

* fe (*i,* 2 ,Ai,A 2 ) = 0, *;(ti, t 2 , Ai, A 2 ) = 0, (1.5) 

for which (0, 0, A*) is a solution and d(^ k , * ; )/<9(Ai, A 2 ) at (0, 0, A*) is invertible. Then 
the implicit function theorem gives A in a neighbourhood of A* in K 2 as a function of 
(ti, t 2 ) in a neighbourhood of the origin in K 2 . This is the sheet of general solutions we 
arc seeking. It intersects the "special" sheets when t\ = or t 2 = 0. The present analysis, 
which yields a qualitative local description of the set of all 27r-periodic hydroelastic waves 
near a bifurcation point, is sufficient to show that secondary bifurcations occur. 
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Figure 1: A possible bifurcation diagram in the space (ii,i2,Ai), when A2 is fixed. The 
dashed curves correspond to the two branches of "special" solutions on the planes t\ = 
and £2 — 0, and the solid curve gives the secondary branch of "general", symmetry- 
breaking solutions. 

Remark. A more detailed geometrical description of the bifurcating sheets depends on 
the coefficients in Taylor series arising in the bifurcation equations, and to calculate their 
values can be very complicated. For example, suppose that we want to draw a picture of 
the solution set when one of the parameters, A2 say, is fixed, and assume that (|1.5p has 
the form 

/(A x ) + {At\ + Bht 2 + Ct\) - A 2 = 0, 
ff(A0 + [at\ + I3ht 2 + 7t|) - A 2 = 0, 

where / and g are smooth functions and A, a, B, 0, C, 7 are constants. Suppose also 
that 

Then it is clear from the implicit function theorem that locally 

Ai = A((A - a)t\ + (B- P)t!h + (C - T)*l), 

for some function A with A'(0) 7^ 0. Therefore, for fixed A2 close to A 2 the solution set 
is given locally by the level set of the function 

f(A({A - a)t\ + (B — 0)ht 2 + (C- 7)^)) + At\ + Bt x t 2 + Ctj. 

The complexity of the dependence of this set on A, a, B, (3, C, 7, / and g is evident. 
These quantities in turn depend, in an explicit but highly non-trivial way, on the elastic 
properties of the membrane. □ 



S 



Remark. In the case when max{fc, I}/ min{fc, 1} G Z, analysis is still possible, but the 
details are yet more complicated. □ 



Remark. In the present work solutions will be found for values of parameters that are 
not covered by the maximization argument in [7j. Indeed, a convexity hypothesis of the 
form En(v, fx) > c 2 p = Ai was crucial in the existence proof of [7], whereas both primary 
and secondary bifurcations occur here provided only that Ai ^ E±±, see Lemmas 15.11 and 



2 Mathematical formulation 

Suppose that in the free-boundary problem (|1.3[) . the shape of S is known. Then ip is 
given by the unique solution to (|1.3l d.e.f). Thus the kinetic and potential energies of 
the fluid are determined solely by the shape of S. On the other hand, the elastic and 
gravitational potential energy of the membrane are determined by the positions of the 
material points (2, 0) in the deformed membrane. To deal with this distinction, and 
ultimately to prove in Section [4] that only the shape matters, suppose that the shape of 
<S is given by a parametrization 



where x '■ ^ ~^ K is a diffeomorphism with x( s + 27r) = 2n + x(s). As before, s is 
the steady travelling- wave Lagrangian coordinate, the unknowns are g(r) and x(s), and 
r(s) = g(x(s)), so that S = {r(s) : s € M}. To develop this approach, we recall the class 
of parametrizations in [7] , specially tailored for this problem. 

Let wbea 27r-periodic real- valued function with second derivative locally square- integrable 
on R, and let C denote its Hilbert transform. Then 



is a 27r-periodic curve in the plane. Thus w is the unknown that describes the wave 
shape. The other unknown is the stretch of the reference membrane. To describe it we 
follow [2] by introducing diffeomorphisms k(t) (k = x^ 1 in (|2.ip ) of the interval [0, 2ir] 
such that x(0) = and x(2tt) = 27r. Then if the material point s of the membrane in 
the reference configuration is 
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□ 



S = {q(t) : t G R}, where q{t + 2tt) = (2tt, 0) + q{t). 
Then, as in [7], we seek R for travelling waves in the form 

R(x, t) : = (cot, 0) + g(x(x - ct)), 



(2.1) 



q(w)(t) := (-T -Cw(t), w(t)), tGM, 



(2.2) 




(2.3) 
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where 



n(w)(r) = Vw'(r) 2 + (1+Cw'(r)) 2 . 
Also, since r(s) = q(w)(t), there is a useful formula for the curvature, 

1 „/fi(«0'(7t 



il{w)(T) V U{w)(t) J 



Roughly speaking, given a profile parametrized by (|2.2j) . the diffeomorphisms (|2.3|) de- 
scribe the family of all physical deformations (s,0) i— > r(s) o/ f/ie reference state which 
produce the same profile £>(w) . Note that the parametrization of a curve and the stretch 
of the material in it are independent. It will be convenient later to replace k with the 
new unknown 

£(r) := «'(r) - 1. 



Then £ has zero mean, 



i r 2 " 

sj[ 5(T)dr = °- 



because «(0) = and k(2tt) — 2ir. In [7] it is shown how solutions to this hydroelastic 
wave problem with a heavy membrane corresponds to critical points of a Lagrangian 
which, in terms of the unknown (w, £), is written 

r 2 r 27T a f 27T 

J(w, £):=-£/ wCw'dT-^l w 2 (l+Cw r )dT (2.4) 

p 2tt 



2 Jo 2 Jo 



- gp / + C) dr. 

Jo 

The first term in the right in (|2.4[) is the fluid's kinetic energy in one period, relative to 
the moving frame; the second term, with a minus, is the change in gravitational potential 
energy of the same body of fluid relative to a uniform flow; the third is minus the elastic 
energy of one period of the deformed membrane; the fourth is plus the kinetic energy of 
the membrane; the fifth term is minus the gravitational potential energy of one period 
of the membrane. (For the derivation, see the discussion leading to [7, (2.18)].) 

This will be the starting point for this analysis. Local bifurcation theory will be used 
to give a complete descriptions of all small-amplitude 27r-periodic waves represented by 
critical points of J close to a bifurcation point. 

3 The equations 

Define the notation (E — V&)(f,/() by 

{E - VE-) (i/, n) := E(y, /i) - v E t (u, (i) - ft E 2 (u, l x). 
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Suppose that w and £ are small in an appropriate norm and that J is differentiable at 
(w,£). Then the partial derivative with respect to £ in a direction 77, where [rj] = 0, is 

. t/ * /" 27r f/r, v^/^H n(«;)o-(iD)\ cV/%)\ 2 1 , 

Let £[u>], depending on iu, be the linear operator defined in [2J Section 4.5] by 

„, w > w'u + (1 + Cw') Cu 

c[w]{u) : = o^p ' 

with the property that 

d w n(w)h = n(w)C[w](ti), d w (n(w)a(w))h = -C(C[w](h'))'. (3.1) 
Then the partial derivative of J(w, £) with respect to w in the direction h is 

,-2k 



(It 



d w J{w,^)h = j h |cq Cw' — gw(l + Cw') — gC(ww') — gp(l + £)| dr 

+ £" £3 (2M.2MfM)( C£ H(V))'^ 

Now, suppose that (u>,£) is a critical point of J. Define the projection 

2 



Pit := u — — I u(t) dr, 
27r Jo 

for all 27r-periodic, locally integrable functions u. Then Vu has zero mean on [0,27r]. 
Note that the operator C[w] is independent of the mean of w. A simple calculation 
shows that d w J(w, £) = if and only if 



1 

— (w + wCw')dT + p = Q and d w J(w,£,)Ph = V h. 
27r Jo 

Hence it suffices to study the equation dJo(w,£) = 0, where Jo is defined by 

( 1 f 2?T \ 2 
Jo(w,£) ■■= J(w,£) +gir ^— J wCw' dr + pj , 

over a class of functions satisfying 

I w(r) dr = £(r) dr = 0, 
Jo Jo 

because then dJo(w, £) = implies that dJ(w*,^) — 0, where 



1 /- 27r 

10* := w; / w Cw' dr — p. 

2^ Jo 
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From a calculation similar to that for J, partial derivatives of Jo are given by 

(3.2) 

and 



d w Jo(w, £)h = j h ^CqCw' — gw(l + Cw') — gC(ww') — gp(l + dr 

f 2 * ( 2 n( w ) _ /f2(to) n(to)(7(io)\\ , „. 
+ j / o (c 2 P T ^-^( T ^ ) ^ r ^i))^)£ M (/ i ')^ 

+ 2g ('hJ wCw ' dT + P)J Q hCw'dr. (3.3) 

In [21 Sections 4.1 & 4.2] the membrane density p is zero. However, the calculations 
there are easily extended to take account of the extra terms here which involve p > 0. 
To proceed, we adapt the notation from [2] to the case of positive p. For u with zero 
mean, let 

VI := ( c o + ^ J wCw ' dT + 2 9P) Cw ' - 9 w i l + Cw ') - gC{ww') - gp(l + f). 

With the L 2 -adjoint of the inverse operator £[u;] _1 given by 

(£[iu] _1 )*(ti) = w'u + C((l + Cw')u), 

let 

mo-iCiw]- 1 )*^^ P(VJ )tft). (3.4) 
Then, as in [2], the equations for critical points of Jo can be written as follows: 
„f , „„,/fi(«/) f2(uA<j(w;) \ c 2 p /Q,(w)\ 2 } 



^)*} = 0. (3.5b) 



2 ^H' 



1!' 

Note that (w, £) = (0, 0) solves (I3.5P for all values of c, c , g, p. 

The free-boundary problem in Section 11.11 is for symmetric waves, so we can simplify 
matters by studying the bifurcation problem in spaces of even functions. This is what 
we do subsequently. Note that, if w and £ are even functions, then Cw' , Q(w), <t(w) and 
1 + £ are also even. 
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Lemma 3.1. Suppose that (w,£) are even functions such that 

dJ (u>,O(M)=0 (3.6) 
for all even {h,rf) with sufficient regularity. Then (13. 6ft holds for all (h,T)). 

Proof. Since every 27r-periodic function is the sum of even and odd functions, by the 
hypothesis it suffices to observe that (|3.6p holds for (h,rf) odd when w, £ are even. This 
follows by (|3.2[) and (|3.3|) , since C[w](h') is odd for odd h. □ 

4 A further simplification 

In this section we use the implicit function theorem to show that, for solutions of (|3.5a[) . 
£ is a function of (w, Ai) near w = £ = 0. This means that the stretch variable £ can be 
eliminated and the problem reduced to one for the unknown shape which is given by w. 
For k £ N, let Hq denote the space of real- valued, 2-7r-periodic, even, zero-mean functions, 
with fcth weak derivative locally square-integrable. For r > 0, let B r (X) denote the open 
ball of radius r centred at the origin in a Banach space X. Fix r > such that 

\<n(w)<2, \a(w)\<l, |£| <I 

for all (w,£) £ B r (H^) x B r (H%). Then a map M: B r (H%) x B r (H%) x (0, +oo) -> i^ 1 
may be defined by 

, , n ^f,„ „^,/fKw) n(w)a(w)\ Ai /f2(u>)\ 2 ~l 

because £7 is smooth, and £l(w)/(l + £) and fl(w)a(w)/(l + £) are bounded functions. 
Thus, the Euler equation (|3.5a[) may be written as 

M(w, £, Ai) = 0. 

Moreover, M(0, 0, Ai) = for all Ai. We note that 

M e C ,oc (B r {H^) x B r {Hl) x (0,+oo), Hi), 

and, when (w, £) = (0,0), 

n(0) = l, cr(0) = 0, dO(0)ft = C/i', d(Ocr)(0)ft = /i". (4.1) 

Since £[0] = (C[0]- r )* = C, it follows from F^and (33} that 

4M(0,0,Ai)/i, = -(£ u - Xi)Ch' + gph, (4.2) 
d e M(Q,0,Ai)7j = (£;ii-Ai)ry, (4.3) 
d Al M(0,0,Ai) = 0, (4.4) 

where for convenience we have written En instead of .En (1,0). 
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Lemma 4.1. When Ai ^ En, there exist a neighbourhood W of (0, Ai) in Hq x (0, +oo) 
and a map £ G C°°(>V, B r (H„)) such that 

M(w,£(w, Ai), Ai) = 

/or a/Z («;,Ai) G W, and, i/ M(w, £, Ai) = (u>,Ai) £ W, i/ien £ = ?(w,Ai). 
Moreover, for all Ai 7^ Bn, 

!(0,A 1 )=0, (4.5) 

c«(0, Ax) ft = Cti + - 9f } - h. (4.6) 



Proof. To obtain the existence of £, apply the implicit function theorem using (|4.3[) . The 
required formulae then follow by the chain rule and (|4.2p and (|4.4[) . □ 

Because of this, all solutions (w, £, A) of (|3.5p with (w, Ai) 6 W are solutions of (|3.5b() of 
the form (10, £(w, Ai), A). Let to(u>, A) denote mo (see (|3.4[) 1 when £ = £(u;,Ai), and let 

eifiw.Ai) := J5i( =^ , — _ v ; ), i = l,2. 

Then on the set V := {{no, A) : (w, Ai) G VV, A2 > 0}, define the function 

Ai^M 2 



F(w,X) := Ve 2 (w ) Xi) +C\ / r(m(w,X) + n{w)e 1 (w,X 1 ) - ) dt\, 

wo v l + £(w,\i)' > 

so that the system (|3.5[) . for (id, Ai) G W and Ai ^ En, becomes 

F(w,X) = 0, (w,X)eV. (4.7) 



5 The linearized equation 

Recall that F(0, A) = for all A with A x ^ E n and that F G C°°(V, Hi). Wc now 
calculate its partial derivatives. When w = 0, (|4.1[) holds, £ = by (|4.5|) and, as a 
consequence, V(VIo) = 0. Hence, by (|4.6[) . 

^to(0, A) ft = -(A 2 + gp) h + ( ^ P)2 Ai -9) C (J q h(t)dt), 

and, by H2, 

d w F{0,X)h= E 22 h" + Xxh- X 2 c(J h(t)dt) (5.1) 
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where, as above, En — Ea(l, 0), i — 1, 2. 

Now suppose that h G Hq\ {0} and that d w F(0, \)h = 0. Then h e H%, and differenti- 
ating the equality d w F(0, A) h — twice with respect to r yields 

E 22 h"" + X 1 h" - A 2 Ch! + (g + ^ 9P) 1 ) h = 0. (5.2) 

V Ax — Exi ' 

It is easy to see that (|5.2[) has a non-constant even solution h if and only if 

E 22 k A - Aifc 2 - A 2 fc + g + (gp) F = 0, (5.3) 

Ai — i/n 

for some positive integer fc. Since -££22(1,0) is assumed to be positive, for every fixed 
Ai, A2 > 0, (|5.3p possesses at most two positive integer solutions. (This follows by noting 
that the graph of the quartic curve x E 22 x a — \\x 2 on the half-plane {x > 1} intersects 
any straight line with slope A2 at most twice.) 



5.1 Non-trivial kernel 

Let g, p, En,E 22 > be fixed. Then Ai, A 2 > 0, Ai — En ^ and (|5.3f) holds for some 
integer k > 1 if and only if 

< \ 2 k = E 22 k 4 - Axfc 2 + g + ^ 9P)2 = ^ Xl \ and Ax > 0, (5.4) 

Ai — E\\ En — Ai 

where 

p k (X) := k 2 X 2 - (E 22 k 4 + Enk 2 +g)X + En(E 22 k A + g) - {gp) 2 . 

The discriminant of pk is 

A(p fc ) = (Enk 2 - E 22 k 4 - gf + (2gpkf > 0, 

and, since pk(En) = —{gp) 2 , its roots 

v ± _ Enk 2 + E 22 k 4 +g± y/Afa) 
Xk - 2k 2 

satisfy 

*k < E U < x+. 

Moreover X^ > if and only if En(E 22 k A + g) > g 2 p 2 . It follows that (|5.4j) holds if and 
only if 

Xi e(0,X fc -)u(£u,X+), A 2 = / fc (Ai), 
where (0, X^) is meant to be empty if X^ < 0, and 

fk(Xi) ■= e 22 fc 3 - \ik + y(g ' {gp)2 



k V Ai - E u 

X7 — ► En and X^ — > +00 as — * +00. Thus, for every Ai 7^ -En, there exists an 
integer k = fe(Ai) such that Ai 6 (0, -XT) U (£11, Xt) for all fc > fc. Therefore, for every 
fc > fc, (|5.3|) holds with A2 = /fc(Ai), and we have proved the following lemma. 
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Lemma 5.1. Let g, p, En and E22 be fixed, positive constants. For every fixed X\ 7^ -En, 
the parameters A2 for which the linearized operator d w F(0, Ai, A2) has a non-trivial kernel 
form a sequence {X^ = /fe(Ai) : k > fc(Ai)}, with 

A 2 — /fe(Ai) — > +00 as k — > 00. 

Thus, for every g, p, En, E 22 > there exists a set A formed by infinitely many curves 
Ak in the parameter quadrant {(Ai, A2) : Ai > 0, A2 > 0}, 

A={jA k , A = {(Ai,A 2 ):A 2 = / fe (Ai)}n{Ai>0, A 2 > 0} 

(see figure^]), such that the kernel of the linearized operator d w F(0, A) is nontrivial if 
and only if A S A. Note that there is no restriction on Ai except that Ai ^ En- 

5.2 Double eigenvalues 

The kernel of d w F(0, A) is two-dimensional if and only if (|5.3[) has two positive integer 
solutions k ^ I, namely the curves Ak and Ai cross at A. Now, k ^ I solve (|5.3p if and 
only if 

A 2 = fefc.i(Ai) := {k + l)(E 22 (k 2 + l 2 )- Ax) and g M (Ai) = 0, (5.5) 

where 

q M pf ) := MX 2 - (Enkl + E 22 kl(k 2 + kl + Z 2 ) - g)X 
+ E 11 E 2 2kl{k 2 + kl + I 2 ) - E n g + [gpf- 
For all kl sufficiently large, the discriminant of qk,i, 

A(q k j) = (Enkl - E 22 kl(k 2 + kl + I 2 ) + gf - 4kl(gp) 2 , 
is positive, and the roots XZi and Xt, of qk,i are both greater than En- 

Since hk.iiX^) < for all kl sufficiently large, there are at most finitely many solutions 
A of ([S3)) with Ai = X^ and A 2 > 0. On the other hand, hk,i{XZ{) > for all kl 
sufficiently large, and 

as kl — > +00. Thus, we have proved the following lemma. 

Lemma 5.2. Let g,p,E\\ and E22 be fixed, positive constants. The parameters A for 
which the linearized operator d w F(0, A) has a two-dimensional kernel form a sequence 
AW = (\[ n) ,X { 2 n) ), with 

\ (") \ jy \ («) 1 

Aj \ En, X 2 — > +00 as n — ► 00. 

Remark. Double eigenvalues with Ai < En are possible, provided we assume some 
additional hypotheses on E, namely 

E 2 2{kl) 2 <g< E 22 kl{k 2 + kl + I 2 ) 

and En sufficiently large. In any case, they are at most finitely many. □ 
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6 Lyapunov-Schmidt reduction 



We turn now to study the bifurcation of solutions of (|4.7D . Recall that throughout we 
are dealing with 27r-periodic functions w of zero mean. Suppose that A* = (A*, A2) G A, 
XI ^ En. Then the kernel, 

V := Kevd w F(0, A*) C H%, 

of the linearized operator, is a subspace of dimension 1 or 2, depending on the number 
of integer solutions of (|5 ,3|) , and the range 

R := Range ^^(0, A*) C Hi 

is orthogonal to V with respect to the L 2 (Q, 2ir) scalar product, namely 

Hl = V®R, Hi = V®(RDH^). 

In fact, it is evident from (|5.ip that d w F(0, A) is a diagonal operator with respect to the 
basis of even 27r-periodic functions {cos(jr) : j = 1, 2, . . .}, for all A. 

Following the classical Lyapunov-Schmidt decomposition, we write 

w = v + y, veV, yeRnHQ, 

and denote Hy,!^ the projection onto V and R respectively. The equation F(w, A) = 
is then equivalent to the system 

Ily F(v + y, A) = (bifurcation equation), 

(6.1) 

IIr F(v + y, A) = (auxiliary equation). 

Lemma 6.1 (Auxiliary equation). There is a neighbourhood U of (0, A*) in V x M 2 , a 
neighbourhood U of in RCi Hq and a function y S C°°(14, U) such that 

tt R F{v + y(v,\),X) =0 

for all (v, A) G U, and, if lip, F(v + y, A) = with y G U and (v,X) G U, then y = y(v, A). 
Moreover, for all (0, A) G U, 

y(0,A) = 0, d„y(0,A)=0, d Ai y(0,A) = 0, i = 1,2, (6.2) 

and t/iere exists a constant C > suc/i t/ia£ 

II^.AJIIh, <C||t;|| 2 ff3 

for all (v, A) G U. 

Proof. Apply the implicit function theorem, and note that, for all A, d w F(0, A) is diagonal 
in the basis {cos jr}, j = 1,2,.. ., therefore n# d w F(0, X)v = for all v G V, for all A. □ 

In this way, the bifurcation problem for the equation (|4.7|) has been reduced to 

Il v F(v + y(v,X),X) =0, (6.3) 

with (v, A) G U C V x M 2 . 
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7 Sheets bifurcating from a simple eigenvalue 



Here we study the elementary case in which equation (|5.3[) . and hence the linearized 
operator (|5.ip . has a 1-dimensional kernel. From the discussions in Section [5] this is the 
case for parameter values (AJ, Aj) on the union of countably many curves with countably 
many points removed. 

Hence suppose that there exists a unique integer k > 1 that satisfies (|5.3[) for A* = 
(Ai,A£). Then 

V := Kerd w F(0, A*) = {tcos(kr): t G M} 
and, by (|6.3p , the system (|6.1[) is equivalent to the problem 

$0,A) := ny F(tcos(kr) + y(tcos(kr), A), A) =0. (7.1) 

$ is a smooth real valued map of three real variables (t, Ai, A2), defined on a neighbour- 
hood of (0, A*, A|). Since F(0, A) = for all A, from Lemma [BTTI it follows that 

$(0, A) = for all (0, A) G U. 

Also, using (|6.2p and the orthogonality of V and R, 

a t $(0,A*) = 0. (7.2) 

To find solutions of (|7. 1[) with i ^ we invoke the implicit function theorem in the usual 
analytic approach to bifurcation problems. 

Theorem 7.1. Suppose that there exists a unique integer k > 1 that satisfies (|5.3p for 
X* = (A*, A2), with XI ^ -En- TTien there exist neighbourhoods IA\ of (0, AJ) m K 2 and 
C^i 0/ A2 in K, and a map A2 G C°°(Z/l, t^i); wiift A2(0, A*) = A2, such that 

$(i,Ai, A 2 (*,Ai)) = /or a/Z (i,Ai) eWi, 

cmd, i/ $(£, Ai, A2) = 0, with (t, Ai) G Wx> * and A 2 G U±, then A2 = A2(i, Ai). As a 
consequence, 

F(w(t, Ai),Ai,A a (*,Ai)) = 0, 

w/iere 

iu(t, Ai) := tcos(fcr) + y(t cos(fcr), Ai, A2(t, Ai)) = t cos(fcr) + 0{t 2 ). 
Proof. First, we prove that 

a t 2 A2 $(0,A*)^0. (7.3) 

By (E2D, 

<9 t $(0, A) = d w F{0, A) (1 + d v y(0, A)) cos(fcr) = H v d w F(0, A) cos(fcr), 

and 

9 t 2 A2 $(0,A*)=n v d 2 o>A2 F(0,A*) cos(fcr). 
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By & 



d™F(0, A) cos(fcr) = ^ - fc^ 22 + A i + y - (5 + A J fc2 ) cos(fcT), 

whence 

n v < A2 F(0,A*)cos(fcr) = i >0, 

and (|7.3p is proved. Since 

$(t,A)= / (9f$)(x,A)da; = i / (<9 t $)(zi, A) dz, 
Jo Jo 

it follows that <I>(i, A) = 0, with t ^ 0, if and only if A) = 0, where 

<p(t,\) := [ (d t $)(xt,X)dx. 



From the smoothness of $ it follows that (p is also smooth. By (|7.2p . ^(0, A*) = 0. 
Moreover, d\ 2 ip(0, A*) ^ by (|7.3p . The result now follows from the implicit function 
theorem. □ 

Remark. Since 

H V d 2 wM F(0, A*) cos(fcr) = 1 + (^-^—^ > 0, 
the role of Ai and A2 can be swapped. □ 



8 Bifurcation from a double eigenvalue 

We have observed that for any (Ai, A2) there are at most two positive integer solutions, 
k, I, of (|5.3[) , and this happens only if fk(Xi) = fi(Xi) = A 2 . Suppose that there are 
indeed two such solutions, k and I, with 



max{fc, 1} 
min{fc, 1} 



Z. (8.1) 



Let Zk be the closure of span {cos(j kr) : j £ N} in L 2 (0, 2ir), and similarly for Z\. Now 
note that if one seeks waves with minimal period 2ir/k or 2tt/1, the original bifurcation 
problem (|4.7[) may be specialized to a problem on Zk or Zk and the reduced problem 
(|6.3p is similarly restricted to Zk or Zi . In each of these restricted settings separately, 
only one solution, k or I, of (|5 .3[) is relevant, and there is a simple eigenvalue from which 
a curve of solutions in Z k fl Hq or Z[ n #q bifurcates, exactly as in the preceding section. 
However, we will now show that other solutions that are neither in Zk nor Z\ bifurcate 
at A* when k and I are solutions of (|5.3[) with A = A* and (|8.ip holds. 
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In this case the kernel of the linearized problem is two-dimensional, 

V := Ker d w F(0, A*) = {h cos(fcr) + t 2 cos(/r) : (ti,t 2 ) S M 2 }, 
and the bifurcation problem (|6.3p is 

*(ti 1 i a ,A) = ) A = (Ai,A 2 ), (8.2) 

where 

$(ti, i 2 , A) := + A), A), t> = ti cos(/cr) + t 2 cos(Zr). 

Let $fe cos(fcr) := IIfe$ and $; cos(Zt) := II; where life and II; denote the projections 
onto span{cos(fcr)} and span {cos(lr)}, respectively. Thus (|8.2|) becomes 

$ fe (ii,t 2 ,A 1 ,A 2 ) =0, 
$,(*i> i 2 , Ar, A 2 ) = 0, 

a system of two equations in four unknowns which is satisfied by (0, 0, Ai, A 2 ) for all A. 
The key to our result is the following observation. 

Suppose that t\ = 0, and v — £ 2 cos(7t), t 2 G K. Then an application of Lemma 16.11 in 
the subspace Z; of 27r/Z-pcriodic functions yields that y(v, A) G Zi n R, because of the 
local uniqueness in the implicit function theorem. Hence v + y(v,\) is 27r/7-periodic, 
therefore F(v + y(v, A), A) is also 27r//-periodic. As a consequence, 

$fc(0,t 2 ,A)=0 far all t a , A. (8.3) 

For the same reason, 

$;(ti,0, A) = for all ti, A. (8.4) 
We now require the non-degeneracy condition 

Remark. It is easily checked that condition (|8.5p is equivalent to the geometrical as- 
sumption that /fc(A*) 7^ /;'(A*). In other words, the curves ^ and Ai are not tangential 
at their intersection point A*. □ 

Theorem 8.1. Suppose that there exist two integers k,l that satisfy (|5.3I) for A* = 
(A|,A 2 ) where XI =/= En, and ()8.ip and (|8.5[1 ZioZd. T/ien i/iere exist neighbourhoods 
Ui of the origin and C/ 2 of A* m R 2 , and functions A(ti,t 2 ) = (Ai(ti, t 2 ), A 2 (ti, t 2 )), 
Xe C 0O (U2,U 2 ), with A(0, 0) = A*, sttc/i i/iai 

$(ti,*a,A(ti s t a )) = /or oH (ti,t 2 ) G £Y 2 , 

and, i/$(ti,t 2 ,A) = iwt/i (ti,t 2 ) G W 2 \ {(0,0)} and A G U 2 , then X = A(ti,t 2 ). .As a 
consegwence, 

F(^(ti,t 2 ),A(ti,ta)) =0, 

where 

ttf(ti,ta) := ti cos(fcr) + t 2 cos(Zr) + y(ti cos(fcr) + t 2 cos(Zr), A(ti,t 2 )J 
= ti cos(fcr) + t 2 cos(Zr) + 0(i 2 + f 2 ,). 
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Proof. Let # := (* fc ,*0) 

#fc(*i,*2, A) := ^ (a tl $fc)(xfi,t 2 ,A)dx, #i(*i,*2,A) := / (<9 t2 $ fc )(*i> xt 2 , A) dx. 
VP/t and are smooth by (|8.3p and (|8.4p . Moreover, since 



* fc (0, 0, A) = n fc d w F(0, A) cos(fcr) (8.6) 

= -E 22 k 2 +X 1 + ^L -fg- 



^ .hi ( 

^2 ' 



Ai - En) k 



and, analogously, 

*«(0,0,A) = —E->>f- ~ A i + -^ - (//-- x )-• 1^.7) 

it follows that 



2 ' A)?- 



#(0,0, A*) = 0. 

To apply the implicit function theorem to # at the point (0,0, A*), it is sufficient to 
prove that the 2x2 matrix representing the linear map d\^>(Q, 0, A*) is invertible. Now, 
differentiating (18. 6p and (18. 7|) with respect to Ai and A2, 

I 1 \ 1 
~kV 



det 

which is nonzero by (18. 5|) . □ 



Remark. By the definition of #, for (iijia) 6 W2, with ii ^ and £2 7^ 0, Theorem 18.11 
gives solutions of problem (|8.2[) which do not belong to nor Zi, as it had been stated 
above. □ 
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